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Introduction 

The harmonic series 1 + 1/2 + 1/3 + 1/4+.. . is infinitely large, whereas the series 
of squared terms 1 + 1/4 + 1/9 + 1/16 + ... sums to a finite value. Indeed, the 
divergence of the harmonic series is easily verified by grouping its terms in packages 
whose lengths are increasing powers of two: 

1 /l 1\ /l 1 1 1\ 12 4 

1+ 2 + U + i/ + U + 6 + 7 + 8j + --' > 1+ 2 + 4 + 8 + - 

and the boundedness of the squared series is equally simple to check: 

111 111 

1+ 4 + 9 + l6 + -" < 1+ l-2 + 2^3 + 3^ + --- 

Alternatively, one can of course also establish these well known facts by comparing 
with integrals: 

Ai r" dx , A i r" dx i 

n— 1 u n— 1 1 

The problem of computing the exact value of 1/n 2 is so classical (dating 
back to 1644) that it even has a name of its own: the Basel problem. The city of 
Basel was the hometown both of the famous brothers Jakob and Johann Bernoulli, 
who made serious but unsuccessful attempts to solve the problem, and to the prodi- 
gious mathematician Leonhard Euler, who found the first solution in 1735, see [3]. 
Since then a great many different ways of evaluating the sum have been discovered. 
Fourteen proofs are collected in [2] , and two of them have been in included in pQ . 

In this note we shall present a new simple geometric method for finding the 
exact value of Y^nLi V^ 2 using essentially just basic trigonometry: the sine rule for 
triangles; combined with some elementary undergraduate analysis: area-preserving 
maps, see for instance [3]. 
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Bipolar coordinates and the (co)sine rule 

Consider an arbitrary triangle with base of length 1, such as the one in Fig. 1 below. 
The position of the top of the triangle can then be determined either by the two 
remaining side lengths A and B or, equivalently, by the two interior base angles a 
and (3. One can view these data as representing two alternative sets of coordinates 
for the point at the top of the triangle. Thinking of the two base vertices as a pair 
of reference points (much like your own two eyes) it is natural to refer to (A, B) 
and (a, (3) as bipolar coordinates, radial and angular respectively. 
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FIGURE 1. The side lengths and interior angles in a triangle. 

Let us check what are the possible values of the two sets (A, B) and (a, (3) 
of bipolar coordinates. For the radial coordinates (A, B) the triangle inequality 
amounts to the obvious restriction that each side length must be less than the sum 
of the remaining two. That is, one has the following three inequalities: 

1<A + B, A<l + B, B<l + A. 

These are in fact the only conditions imposed on the radial bipolar coordinates, so 
the collection of all possible (A, B) is given by the infinite polygon S depicted on 
the left in Fig. 2. For the angular coordinates {a, (3), which we always measure in 
radians, the corresponding restriction comes from the fact that the sum of all three 
angles in a triangle is equal to it. The possible values of the angular coordinates 
(a, (3) are therefore given by the half square T shown on the right of Fig. 2. 




Figure 2. The possible values for radial and angular bipolar coordinates. 
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To solve a triangle means to find all its side lengths and all its interior angles, 
with only some of these data being known beforehand. This is precisely what one 
needs to do in order to pass from one set of bipolar coordinates to another. The 
classical tool for this is the sine rule 



sin a 
~A~ 



sin (3 sin(7r — a — (3) 



B 



1 



sin(a + 0) , 



(1) 



which permits one to express the side lengths (A, B) in terms of the angles (a,/?), 
and the cosine rule 



A = 1 + B — 2B cos a , B 2 = 1 + A 2 - 2 A cos (3 , 



(2) 



which conversely expresses (a, (3) in terms of (A,B). More precisely, the sine rule 
provides us with the explicit bijective map F : T — > S given by 



F(a,P) 



sin a 



sin/? 



sin(a + (3) ' sin(a + (3) 
whereas the cosine rule allows us to write down the inverse map 



F-\A,B) = 



arccos 



I- A 2 



2B 



arccos 



1 + A 2 - B' 
2A 




Figure 3. The set of values of the logarithmic side lengths: an amoeba. 

It will now be useful to make a logarithmic change of scale for the side lengths 
in Fig. 1, and thus to replace the radial bipolar coordinates (A, B) by the new 
coordinates (x, y) given by x = — log A, y = — logB. The image of the set S under 
this change of coordinates, that is, the result of drawing it on a log-log graphpaper, 
is the set U shown in Fig. 3. This set U, with its "tentacles" going off to infinity, 
is an example of what is nowadays known as an amoeba in complex geometry, see 
[5] and [BJ. In the first quadrant its boundary is given by the curve e~ x + e~ y = 1, 
which can also be represented as y = — log(l — e~ x ). 
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The advantage of having passed to logarithmic coordinates becomes apparent 
when we consider the composed bijective map G : T — > U given by 



G{a,(3) = 



log 



sin(a + (3) 



sin a 



log 



sin(a + /?) 
sin/3 



Theorem: The Jacobian deteminant of the map G is identically equal to 1, that 
is, G is an area-preserving map. 



Proof: Componentwise differentiation of G yields the Jacobian determinant 

= — cot(a + /3)(cot a + cot 0) + cot a cot /?, 



cot(a + j3) — cot a cot(a + (3) 
cot (a + (3) cot(a + f3) — cot (3 



and this expression is indeed identically equal to 1, in view of the addition formula 
cot(o! + (3) = (cot a cot (3 — l)/(cot a + cot (3). 

Since G is a bijection, that is, a one-to-one map, the theorem has the following 
immediate and remarkable consequence. 



Corollary: 



Area([/) = Arca(T) = tt 2 /2. 



The threefold shape of the amoeba U suggests the one should let its asymptotes 
divide it into three parts Uq, U\ and U2, as indicated on the left in Fig. 4. The 
points (x, y) on either of the asymptotes correspond precisely to the side lengths 
(A,B) ~ (e~ x 1 e~ v ) of the isosceles triangles, for which the interior angles (a,/3) 
satisfy one of the conditions a = (3, a = ir — a — (3 or f3 — n — a — [3, and these 
are just the equations of the three medians of the triangle T . (In particular, the 
origin (x,y) = (0,0) represents the equilateral triangle, with a = (3 = n/3.) In 
other words, the map G sends the medians of T to the asymptotes of U, and one 
obtains a corresponding polygonal subdivision of T into three parts T , T\ and T2, 
shown on the right in Fig. 4. 




Theorem: The area-preserving bijective map (x,y) (—y,x — y) permutes the 
amoeba subsets cyclically: Uq 1— » U\ 1— > U2 1— > Uq- 
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Proof: To check that Uo i— ► £7i one can simply write down the three defining 
inequalities for J7o in the new coordinates (£, 77) = (— y, a; — y): 

x>0^^<?7, y>0 ^ £<0, 1 < e~ x + e~ y <!=> < e~ J ' + 1 , 

and observe that the new inequalities are precisely the ones that define U\. The 
rest of the theorem is proved analogously. 

It follows immediately from this theorem that Uq, U\ and £/ 2 all have the same 
area, which must then be equal to 7r 2 /6, the total area of U being 7r 2 /2. So, either 
using again that G is area-preserving and that G{Tq) = Uo, or directly calculating 
the area of the simple polygon To, one obtains the following conclusion. 

Corollary: Area(£/ ) = Area(T ) = 7r 2 /6. 

Notice that the points in T and Uq correspond to triangles in which the longest 
side is the base of length 1 . 



Spreading and piling of squares 

Geometrically, the Basel problem amounts to computing the total area of a collec- 
tion of squares with decreasing side lengths 1, 1/2, 1/3, 1/4, and in order to 
acheive this, it will be advantageous to spread out each such square by means of an 
exponential function. More precisely, one observes that a square with side length 
1/n has the same area as the part in the first quadrant lying under the exponential 
curve y = e~ nx /n, see Fig. 5. 




Figure 5. The two areas are the same. 
Indeed, the area under the curve is readily calculated by means of the integral 



1 r 

n Jo 



' dx = 



1 



1 



(Incidentally, this elementary fact can also be deduced by observing that the area- 
preserving bijection (x,y) ((1 — e~ nx )/n,ye nx ) maps the area under the curve 
y = e~ nx jn to the square < x,y < 1/n.) 

The next step will be to pile these spread out squares on top of each other. In 
order to understand what such a process will yield, let us take a new look at the 
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harmonic series. Even though the harmonic series itself is not convergent, it can be 
made to converge by introducing powers in the numerators. 

Theorem: For any positive real number t < 1 the power series t + t 2 /2 + t 3 /3 + 



t 4 /4 + 



convergent and its sum is equal to — log(l — t). 



Proof: Differentiation gives 

2 + ~3 + ' ' ' + iV , 



D(t + 



1 + t + t 2 + ...+t N ~ 1 



l-t 
l-t 



N 



which means that 



t 2 t 3 t N 



1 



io I - ' Jo 

Here the first integral on the right has the desired value — log(l— t) while the second 
integral is less than 



dT 



1 



* T N dT 



1 



t 



N 



dT 



which tends to zero as N 



io 

CO. 



f 



-t N log(l - t) . 



Writing t = e~ x , so that t n = e~ nx , we thus see that for every x > the series 
^2 e~ nx jn has the sum — log(l — e~ x ). But piling the spread out squares on top of 
each other, as indicated on the right in Fig. 6, precisely amounts to summing the 
functions e~ nx /n, so from the theorem one deduces the following fact. 

Corollary: The infinite pile of spread out squares exactly covers the set Uo defined 
by the inequalities x>0,y>0,y< — log(l — e~ x ). 




Figure 6. Building two infinite piles. Equally shaded areas are the same. 

Since the convergence is not uniform in a; a little argument may be in order 
here: One easily verifies that for any compact subset K of Uq there is some finite 
pile < y < e~ x + e~ 2x /2 + . . . + e~ Nx /N which contains K. On the other hand, 
it is clear that every finite pile is contained in IV 

Conclusion: The series 1 + 1/4+1/9 + ... of squared reciprocals of positive 
Mergers has the exact sum ir 2 /6. 
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In order to emphasize the conciseness of the solution to the Basel problem that 
we have just presented, let us recapitulate it again in formulas: 

— = / dx = — l logfl - e~ x ) dx = I dxdy = / dad/3 = — . 

ti n2 ~i Jo n Jo J Uo J Ta 6 

Here the notations Uq and T again refer to the sets appearing earlier, see Fig. 4. 



Epilogue: the complex logarithm explains the mystery 

The most surprising part of the arguments in the preceding sections is probably the 
fact that the map G, connecting the logarithmic side lengths (x, y) and the angles 
(a, (3), turned out to be area-preserving. Why on earth should the Jacobian of G 
be identically equal to 1 ? 

In order to shed some light on this enigma, we shall take a look at logarithms 
of complex numbers. Vaguely put, it is their marvelous property 

logz = log |z| -Margz, 

combining (the logarithm of) the modulus \z\ with the argument argz, that explains 
the presence of a simple relation between the logarithmic side lengths (x,y) and 
the angles (a, (3). 

Let us now be more precise, and consider the linear equation 1 + z + w = 0, with 
the unknowns z and w being complex numbers. Depicting 1, z and w as vectors 
in the plane, one can interpret the fact that they sum to zero as saying the three 
vectors should form a closed triangle, quite similar to the one in Fig. 1. 
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Figure 7. The triangle formed by the vectors 1, z and w. 

We now get 1 + z + w = 1 + Ae l ^~^ +Be l ^ a+ ^ = 0, by choosing the arguments 
argz and argu> to lie between and 27r, and taking real and imaginary parts of 
this complex equation, we see that it amounts to the two real equations 

A sin (3 — B sin a =0, (3) 

1 - A cos /3 - B cos a = . (4) 

It is an amusing exercise to verify that this system of equations is in fact equivalent 
to the sine and cosine rules (1) and (2). For instance, to deduce the cosine rule 
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from (3) and (4), one can first re- write (4) as A cos = 1 — B cos a. Squaring both 
sides, and using (3) to replace the squared left hand side A 2 cos 2 (3 = A 2 (l — sin 2 (3) 
by A 2 — B 2 sin 2 a, one obtains the cosine rule A 2 = 1 + B 2 — 2Bcosa. 

Letting the arguments argz and argw be more general, it is natural to "lift" 
the complex line l + z + w = 0to the exponential complex curve 

X = {(s,t) EC 2 : 1 + e s + e* = } . 

Denote by P and Q respectively the restrictions to X of the two linear projections 

(s,t) \ — >(Res,Rei) and (s,t) \ — ► (1ms, Imt) 

on the real and imaginary parts, and let U = P(X), T = Q(X) be the corresponding 
images in R 2 . Observe that Q is invertiblc on the interior of T, so we can introduce 
a map G = PoQ^ 1 going from imaginary parts to real parts. This map G is closely 
related to our previous map G. In fact, one has the identity 



G = 



-1 
-1 



oGo 



-1 

1 



Here the two linear maps (represented by their matrices) are both area-preserving 
(since the determinants equal 1), so what remains to be explained is why G has its 
Jacobian identically equal to 1. 

Fix any point pel and let T p denote the (real two-dimensional) tangent space 
to X at p. Since P and Q are (restrictions to X of) linear maps, their derivatives 
P'(p) and Q'(p) coincide with (restrictions to T p of) the maps themselves. The 
Jacobian of G is therefore equal to the determinant of P o Q^ 1 , where we now use 
the same notation P, Q to mean the real and imaginary projections T p — ► R 2 . 

Notice next that T p is in fact a (one-dimensional) complex subspace of C 2 , 
so coordinate -wise multiplication by i defines a linear map J on T p . Since J just 
amounts to a rotation (by an angle 7r/2) around the origin it is clear that J is 
area-preserving and thus has determinant 1. Moreover, one has J = Q^ 1 o P, so 
P o Q^ 1 = Q o (Q^ 1 o P) o Q^ 1 also has determinant 1, which means that G has 
Jacobian equal to 1. 
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